Critical Frontier of the Triangular Ising Antiferromagnet in a Field 
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We study the critical line of the triangular Ising antiferromagnet in an external magnetic field by 
means of a finite-size analysis of results obtained by transfer-matrix and Monte Carlo techniques. We 
compare the shape of the critical line with predictions of two different theoretical scenarios. Both 
scenarios, while plausible, involve assumptions. The first scenario is based on the generalization 
of the model to a vertex model, and the assumption that the exact analytic form of the critical 
manifold of this vertex model is determined by the zeroes of an 0(2) gauge-invariant polynomial in 
the vertex weights. However, it is not possible to fit the coefficients of such polynomials of orders 
up to 10, such as to reproduce the numerical data for the critical points. The second theoretical 
prediction is based on the assumption that a renormalization mapping exists of the Ising model on 
the Coulomb gas, and analysis of the resulting renormalization equations. It leads to a shape of the 
critical line that is inconsistent with the first prediction, but consistent with the numerical data. 

PACS numbers: 05.50.-|-q, 64.60.Ak, 64.60.Cn, 64.60. Fr 



I. INTRODUCTION 

The triangular Ising model with equal nearest- 
neighbor coupling _ftr in a magnetic field has the reduced 
Hamiltonian 



Sk 



(1) 



where Si — ±1, and indicates summation over all 

pairs of nearest-neighbor sites. According to the exact 
solution by Houtappel H of the triangular Ising model 
in the absence of a magnetic field, the antiferromagnetic 
model has no phase transition at nonzero temperatures. 
The ground state is characterized by the condition that 
every elementary triangle contains spins of different signs. 
This constraint still leaves a considerable degeneracy, to 
such an extent that the zero temperature antiferromag- 
net has a nonzero entropy. The ground state appears 
to have interesting properties. It is a critical state as 
shown by exact calculations of the spin-spin correla- 
tion function which appears to decay as a power-law of 
the distance. A nonzero temperature T > destroys the 
critical state: the correlations then decay exponentially. 
However, for sufficiently low T, a sufficiently strong field 
H > induces a phase transition to a long-range or- 
dered state, where the minus spins condense on one of 
the three sublattices. As noted by Alexander the 
threefold symmetry of the ordered phase indicates that 
the transition belongs to the three-state Potts universal- 
ity class. The nature of the transition was confirmed b 
Kinzel and Schick |Jj , using phenomenological scaling 
and numerical transfer-matrix calculations; see also Noh 
and Kim 6] and Tamashiro and Salinas 

The critical line covers an infinite range oi K < and 
H. A preview of our numerical data is given in Fig. ^ 
Since the phase diagram is symmetric in H, we restrict it 
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FIG. 1: Numerical results for the {H,e^^) phase diagram. 
The circles (o) denote the data points for Ising temperatures 
T > and the box (■) the so-called KT point at T = 0. 



to H >0. For K — oo, H ^ oo while 6K + H remains 
finite, the model maps 8] onto Baxter's hard-hexagon 
lattice gas of which the critical exponents are exactly 
known, and they do indeed fit the three-state Potts uni- 
versality class. The asymptotic form of the critical line 
in this lattice-gas (LG) limit is 
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where Cc is the exact critical fugacity calculated by Bax- 
ter 0. 

The critical line also extends to K ^ — oo at small 
fields \H\. The behavior of the critical line in this limit 
has attracted attention because of the above peculiar 
ground-state properties, and the associated analytical 
and computational difficulties. It has been conjectured 
Q that the critical line comes in vertically in the ver- 
sus ^ diagram. In other words, when the Ising tem- 
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perature goes to zero, also the reduced critical field He 
(which includes a factor l/k^T) was supposed to go to 
zero. However, Nienhuis et al. |lfl| | provided evidence 
that He instead approaches a nonzero constant when 
K ^ —oo. This result is based on an exact mapping 
of the zero-ternperature Ising model on a solid-on-solid 
(SOS) model jlll |. Using renormalization arguments, 
Nienhuis et al. obtained several critical exponents associ- 
ated with physical fields. It was found that the reduced 
magnetic field is irrelevant: it does not immediately de- 
stroy the critical state a,t K = — oo. 

This renormalization analysis is not rigorous but still 
convincing. Several of its predictions agree with exact 
calculations ^3 aA. H — 0. The renormalization picture 
has been extended to include a nonzero field H and Ising 
temperature T, as well as next-nearest-neighbor interac- 
tions p!I| • It predicts that for T = the model undergoes 
an infinite-order transition to a long-range ordered phase 
at a finite value of the field H. In the SOS language this 
is a roughening transition, in the universality class of the 
Kosterlitz-Thouless (KT) transition ^1. The character 
of this transition was confirmed [T^ Il4j and located at 
Hkt — 0.266(10) by means of transfer-matrix calcula- 
tions and phenomenological renormalization [T^ . The 
associated finite-size-scaling analysis is problematic be- 
cause of slow convergence due to logarithmic corrections 
at the KT transition point. Such corrections are possi- 
bly a reason why an analysis by de Queiroz et al. [l5l |. 
without such corrections, yielded a result that is not fully 
consistent, namely Hkt = 0.211(7) (for the correct in- 
terpretation of this result it is essential that the field H 
used in Fig. 1 and Table I of Ref. does not contain a 
factor 1/T [13). 

The estimated critical field Hkt at T = appears 
to be much smaller than estimates obtained at T > 0. 
The question thus arises whether the Potts critical line 
for T > connects to the KT point at T = 0. It is 
noteworthy that the renormalization scenario given in 
Ref. 10], which includes next-nearest-neighbor interac- 
tions Knnn, implies that the line of phase transitions 
limiting the ordered phase in the {Knnn,T) plane does 
not connect to the transition line in the (i^nnn, H) plane. 
Thus one may ask the same question for the (-R'nnn, H) 
and the {H, T) plane. An answer to this question is 
provided by renormalization arguments presented in Sec- 
tion ^H] This approach also predicts the analytical form 
of the Potts critical line for T ^ while H remains finite. 

A different approach to find the shape of the critical 
line of an antiferromagnetic Ising model in a field was for- 
mulated by Wu ^3 '^ho noted that these models can be 
mapped on vertex models, and that these vertex models 
have symmetry properties that impose restrictions on the 
analytic form of their critical manifolds. He also noted 
that the critical manifolds of the exactly solved vertex 
models are determined by the zeroes of homogeneous 
polynomials in the vertex weights that are invariant un- 
der the symmetry group of the model. On the basis of 



the assumption that the latter form of the critical sub- 
space also applies to vertex models that are equivalent 
with antiferromagnetic Ising models in a field M, one 
may thus attempt to solve for the unknown independent 
coefficients of the homogeneous polynomial, the number 
of which is dramatically reduced by symmetry restric- 
tions. In actual applications, the number of equations 
is still not enough to solve all unknown coefficients, and 
additional numerical input is required, for instance from 
a numerical transfer-matrix analysis. 

This approach is more ambitious than the renormaliza- 
tion analysis in the sense that its aim is to describe the 
whole critical manifold. It has been applied to the Ising 
antiferromagnets on the honeycomb lattice ^.nd on 
the square lattice ^2^. The transfer-matrix data, with 
accuracies in the order of 10^^°, could be successfully 
described by such invariant polynomials of relatively low 
order. Nevertheless one may remark that these analyses 
did not provide solid evidence for the exact form of the 
critical line of the Ising antiferromagnet. 

Application of this approach to the triangular Ising 
model leads to some additional complications. First, the 
topology of the phase diagram is less simple, which re- 
lates to the fact that the lattice is not bipartite. Sec- 
ond, the 3-state Potts character of the critical line implies 
that corrections-to-scaling converge less well in compari- 
son with the Ising case, so that it is not feasible to reach 
the same degree of numerical accuracy. 

In this paper we compare the results of both theoretical 
approaches to our numerical data for the triangular Ising 
antiferromagnet. In Sec. ^ we formulate the invariant- 
polynomial scenario and derive an exact restriction on 
the critical line which must hold if this line is analytic in 
the KT point. We explicitly construct invariant polyno- 
mials of arbitrary order in the Ising vertex weights whose 
roots exhibit this behavior. A summary of the Coulomb 
gas scenario and an analysis of the renormalization-flow 
equations follows in Sec. IIIII The analytic forms of the 
critical lines predicted by these two scenarios appear to 
be mutually inconsistent for T ^ at finite H. In 
Sec. IIVI we outline our transfer-matrix construction and 
present accurate results for the critical points. This sec- 
tion also includes a Monte Carlo analysis of the critical 
amplitudes. An analysis and a discussion of these results 
is given in Sec. Finally, we draw our conclusions in 
Sec. EH 

II. 0(2) INVARIANT POLYNOMIALS IN THE 
VERTEX WEIGHTS 

The mapping of the triangular Ising model on the 64- 
vertex model involves the introduction of bond variables 
bij = or 1 between nearest-neighbor spins and summa- 
tion over the Ising variables ^3 • Since the bond variables 
are independent, there are 2^ = 64 distinct vertices. But 
these turn out to have only 7 distinct weights 

W,,...s,=W, = h''^°''^z'/^ (3) 
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where h = tanh(iJ) and z = tanh(iir). These weights 
are 'symmetric' i.e., depend only on the number of cov- 
ered bonds X^j" ^ij connecting a vertex i to its 6 nearest 
neighbors j. The weights can take imaginary (s =odd, 
K < 0) or real (otherwise) values. The partition func- 
tion of an A^-site system is a homogeneous polynomial of 
order N in the weights: 

^vertex = TT ^si (fc)...^6 (fc) (4) 

Both pertinent indices of neighboring vertices i and j in 
this product are set equal to the bond variable 5^ = 0, 1. 
The summation runs over all possible configurations of 
the bond variables. The partition function of the Ising 
model differs from Eq. Q only by a multiplicative factor 
which is non-singular for finite H and K: 

2cosh(i/)^cosh(if)i^ X Z,,,t,A{Wsih, z)}) (5) 

A crucial property of Eq. is that the summation over 
one of the bond variables 6^ is invariant with respect to 
any 0(2) transformation 

i?(0,£)^f^°%^ (6) 
^ ' ' \^ sm^ ecost^ J ^ ' 

with respect to the indices of the connected vertices. 
Here, det(i?) — e — ±1 distinguishes the S0(2) subgroup 
of proper rotations (e = -1-1) from the improper transfor- 
mations (e = —1) which also include a reflection. Ap- 
plication of this transformation to all bonds connecting 
neighboring vertices (assuming periodic boundary condi- 
tions) leads to a partition sum of the same form but with 
new weights 

W^,,,,,,^R,>,,ie,e)...R,,,,{0,e)Ws,...s, (7) 

where we use the dummy summation convention Si = 
0, 1. This gauge transformation preserves the symmetry 
mentioned under Eq. However, only special 0(2) 

transformations preserve the Ising weight parametriza- 
tion expressed by the right-hand side of Eq. ||2Jl. A triv- 
ial example is the reflection 6 — 0,e — —1 effecting 
Ws{h,z) — > Ws{—h,z). This corresponds to an exter- 
nal field inversion H —H . Below we will first discuss 
another less trivial transformation which also leads to 
weights of the Ising form Q up to a common factor. 
This will have consequences for the asymptotic behavior 
of the critical curve hc{z). In order to investigate the 
assumption that the critical curve is a root of an 0(2) 
invariant homogeneous polynomial in the vertex weights, 
we explicitly construct these polynomials of arbitrary or- 
der. Conclusions about their compatibility with our nu- 
merical data will be drawn in Section Ivl 



A. Dual transformation 

The 0(2) transformation 9 ~ = changes the 
weights as Ws{h, z) z^Wg (± sigii{z)h, 1/z). This con- 
nects (up to a non-singular factor) the Ising weights © 
for physical values of \h\, \z\ ^ 1 to weights of the same 
form but with unphysical values \h\, \z\ ^ 1 correspond- 
ing to complex fields H ± i7r/2, K ± i7r/2. Extending the 
vertex model with weights © to all real values of h, z 
the antiferromagnetic (z < 0) critical curve of this model 
hc{z) has a physical and non-physical branch which are 
connected by the dual transformation 

{h,z)^{h,l/z) (8) 

and which is separated by the self-dual LG and KT 
points at z = —1. This is a rigorous result of seem- 
ingly little use. However, we will now demonstrate that, 
combined with the assumption that the critical curve 
is analytical for z = —1, this severely constrains the 
shape of the antiferromagnetic branches of the critical 
curve for physical values of /i, z. Consider the deriva- 
tives of a branch of the critical line hc{z) at z = — 1, 
assuming it is analytical there. Differentiating the con- 
straint hc{z) — hc{l/z) n times and recursively solving 

for the derivatives hi'^'^ — d'^hc/ dz''\z=-i up to order 
fc = 1, . . . , n one finds that each odd order derivative can 
be expressed as a linear combination of lower even order 
derivatives with integral coefficients. The expansion of 
hc{z) must take the form 

K - hi^~^ + [(1 + zf^ + m(l + z)2™+il + . . . 

2m! 

(9) 

Here to ^ 1 is some integral number i.e., the first non- 
vanishing term is of even order. A notable feature of 
Eq. © is that, if the critical line is analytical at z = — 1, 
it must satisfy dft-c/dz = 0. We now consider the impli- 
cations for the KT and LG asymptotic lines. First, the 
LG asymptotic relation Eq. |2Jl expressed in the variables 
h, z reads: 

which indicates analyticity of the exact form of the LG 
branch of the critical line at z = —1. The first two non- 
vanishing terms in the resulting expansion 

M^) = l-|((H-z)6 + 3(l + z)^ + ...) (11) 

are the same as in the expansion of Eq. (|10|l with m — Q. 
We note that the LG asymptotic curve Eq. (jlUfl is in- 
variant under {h,z) (/i, 1/z) up to all orders. This 
is due to the dual symmetry of the corresponding vertex 
model. Secondly, the approach of the critical curve to the 
KT asymptotic line must also be of the form ^ . How- 
ever, the integer value to ^ 1 in this case is unknown. 
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Expressed in the physical variables H, K, we find a loga- 
rithmic divergence: 



be written in variables s = \ — h? = 1/ cosh^ H and 
M = 1 + z = e^/ coshi^T 



K{H) 



1 

Am 



\n{Hc — Hkt) + const . 



This is the central result of this section. We emphasize 
that it is based on the assumption that the KT branch 
of critical line is analytic at z — —1. Roots of 0(2) 
invariant polynomial equations, which we construct ex- 
plicitly below, all have this property and in the general 
case m = 1. 



B. 0(2) invariant polynomials in the Ising weights 

Now we explicitly construct homogeneous polynomial 
equations in the Ising weights which incorporate all the 
constraints imposed by the 0(2) gauge symmetry (in- 
cluding those discussed above) . Following Perk et al. |23 
we change to the eigenbasis of the SO (2) subgroup of ro- 
tations (e — +1) via the 2x2 matrix Us^ai = ^i—^cnY^'- 



as 



(13) 



Note that = 0, 1 but = ±1. In this basis the trans- 
formation Eq. (UI) takes the simple explicit form (e = ±1): 



A' — pie((Ti-|-...-|-cr6) A 

^(Tl...CT6 ~ ^ ^£Cri...e(T6 



(14) 



For the symmetric vertex model, the 7 components trans- 
form as A'^ = e^^'^A^cr, where a = ai + . . . + = 
0, ±2, ±4, ±6. In this basis the Ising weights are com- 
plex: 



A2k = iz + l) 



2k 



3-fe 



.m=0 



k-1 



m=0 



2k 
2m + 1 



(15) 



where k = 0,1,2,3 and A-2k = ^2fc- Invariant polyno- 
mials which transform with parity ±1 (i.e., I± — > ±/±) 
have the simple form /± i(n^ A^" ± J]^ X"") (or a 
linear combination of these) with exponents which sat- 
isfy aria = 0. A minimal finite set of such exponents 
can be found which generate all other solutions by lin- 
ear combination with integer coefficients. This implies 
that any 0(2) invariant I± can be generated as a polyno- 
mial function of a minimal set of so-called fundamental 
invariants. These 14 polynomials have been constructed 
in pH- The crucial point is to eliminate all dependen- 
cies due to polynomial relations between the fundamen- 
tal invariants (called syzygies) , and further dependencies 
introduced by the parametrization H15|) . To generate in- 
variants of parity -1-1 we need to retain only 4 fundamen- 
tal 0(2) invariants of parity -t-1 which can be compactly 



Iq = ^0 
Ik = A2kA. 



2k 



— ,,6 I ,,6-2/c 



+ u^-^''nkiu)s 



(16) 



where 



Oi(m) = 4(1 -u) 

n2{u) = 4(1 -u)(2u-4)2 (17) 
nalu) = 4(1 -u)(3u-4)2(u-4)2 

The most general invariant polynomial of even order 2e ^ 
2 and parity -1-1 in the Ising weights is 



j=0 



2 



1=1 



(18) 



Cross terms of Ii and I3 have been eliminated using a 
polynomial relation. The (e-l- 1)^ coefficients in Eq. (|18f) 
correspond 1-to-l to an invariant expression in the Ising 
weights of order 2e. To exclude a trivial factorization 
to an even order polynomial, i.e., f2e = -^0/2(6-1)7 we 
require Cg ^ 0. Polynomial invariants of odd order and 
parity -1-1 can be shown to factorize trivially, /2e+i — 
Iof2e, and thus need not be considered further. Finally, 
polynomial invariants of any order and parity —1 can be 
discarded also. We find that for z = —1 such polynomials 
have no other root than h = ±1 i.e., the KT point cannot 
be described. The expansion of f2e in Eq. (fTH|) is not 
well suited to impose restrictions on the coefficients from 
known properties of the critical curve of the triangular 
lattice Ising model. A more convenient but equivalent 
expansion is obtained by replacing Ik — *■ Ik 



r"^ k 

1q , ft 



1, 2, 3 and - 
explicit form 

j=0 



in Eq. (fTH|l . This gives the final 



j 



1=1 



y^<^-3)s3 (19) 



0, s = 1-/ikt 7^ 0) 



Since Ukiu) 42*^-1 in the KT (u 
and LG {s,u — > 0) limit we can easily solve /2e(s,u) 
for the asymptotic relations between s and u in each 
limit. We find that consistent polynomials must satisfy 
Kg — —A^QkI where Cc is given by Eq. (|2J) (LG) and 
Kq = kI = K2 = (KT). This already excludes /2 = 
as a candidate. We thus have (e -I- 1)^ — 5 = 4, 11, 20, 31 
independent coefficients in a polynomial of order 2e — 
4, 6, 8, 10, respectively. The asymptotic value of s in the 
KT limit is determined by skt = 4kq~^/k|. This may be 
used to either extract this value after a fitting procedure 
or as an extra constraint. One can verify from the general 
form Eq. (|19|l that the approach to the KT value indeed 
takes the form Eq. (|12|) with m = 1. This is dictated by 
the dual transformation property Eq. (jS)) combined with 
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the analyticity of the branches determined by f2e = 0. 
Higher integer values for m are also possible but require 
certain coefficients in Eq. H19|l to be strictly zero, which 
is not supported by the numerical data. 



III. RENORMALIZATION ANALYSIS 



A. Mapping on the Gaussian model 

At zero temperature, the three spins of each elemen- 
tary triangle cannot have the same sign. Thus each tri- 
angle has two bonds between antiparallel spins and one 
bond between parallel spins. When all bonds between 
parallel spins are erased, one obtains a lattice tiling with 
rhombi. This tiling can also be interpreted as a stack of 
cubes viewed from the (1,1,1) direction. Thus the zero- 
temperature antiferromagnetic triangular Isingmodel is 
equivalent with a sohd-on-solid (SOS) model [ij. The 
SOS model consists of height variables hi where i denotes 
the lattice site. The height variables assume integer val- 
ues satisfying hi mod 3 = q, where q = 0, 1 or 2 denotes 
the sublattice of site i. Apart from an infinite constant, 
the Hamiltonian becomes 



Ed 



'^|h,-/ij|,l)(l - ^\h,-hj\,2) 



k mod 2,0) 



1) 



(20) 



where we set Kryo — > cxd so that the product (1 — 
^\hi~hj\,i){^ — 5\hi-hj\,2) restricts the height differences 
between nearest- neighbor sites to 1 or 2. The Kronecker 
delta in the second term counts the numbers of + spins. 

The equivalence with the Ising model Q makes it pos- 
sible to express the height-height correlation function [lo| 



^{r) = {[{hr - ha) - {h^ - hr)]'') 



(21) 



in terms of Ising correlations. From the exact results for 
K = -oo, iJ = 0, it follows 10] that 



■ ln(r) 



const 



(22) 



where r is the distance between the correlated sites. This 
result is very useful in the context of a renormalization 
mapping by Nienhuis et al. 'l(f\ on the Gaussian model 
with Hamiltonian 



-H/k^T 



27r 



-h,f- 



2'Khi 
P 



(23) 



where the second summation contains so called spin- wave 
perturbations of the Gaussian model, i.e., a periodic po- 
tential acting on the Gaussian height variables. A term 
with p = 1 originates from the discreteness of the height 



variables in Eq. 1)20(1 . A nonzero magnetic field H favors 
triangles with only one minus spin; in the SOS model 
this leads to an energy alternation between even and odd 
heights. This maps on a spin-wave perturbation with 
p = 2 in the Gaussian model. 

The mapping of Eq. H20|) to Eq. (|23|l is not exact, so 
that the renormalized temperature Tr is in principle un- 
known. However, the height-height correlation function 
of the Gaussian model is known to depend on Tr as 



Comparison with Eq. (|22|l shows that 



Tr = 18 



(24) 



(25) 



for H = 0. Once Tr is known, several quantities of in- 
terest can be calculated for the Gaussian model. These 
quantities include the scaling dimensions associated with 
the so-called spin-wave and vortex perturbations in the 
Gaussian model. Since an elementary excitation of the 
Ising model, i.e., a triangle with three equal Ising spins, 
leads to an SOS height mismatch of 6 units, the Ising 
temperature field ti is represented by the fugacity Vq of 
the q = ±6 vortices. On the basis of the known results 
for the scaling dimensions of Sp and Vq in the Gaussian 
model 0] , a number of properties of the triangular Ising 
model, including a part of the phase diagram extended 
in the direction of Tr, have been derived [lol j. In this 
work we make use of the language of the Coulomb gas 
formulation |2^ to express the relevant scaling dimen- 
sions. The appropriate parameters are the renormalized 
coupling constant ^r, and electric charges e and mag- 
netic charges m. Their relation with the parameters of 
the Gaussian model can be expressed as 



gn = 36/rR 
e — 6/p 

m = q/6 



In this language, the scaling dimensions 
with the activity of charges (e, to) are 



(26) 



associated 



9Rm 



2gR 



(27) 



From Eq. (|25|l we see that g^ — 2 ioi H ~ 0. Since 
the Ising temperature field is associated with magnetic 
charges ±1, the Ising temperature renormalization ex- 
ponent = 2 — Xq,i = 1 is relevant. The system is 
thus disordered for all T > 0. However, the exponent 
yh = 2 — X^ Q = — 1/4 associated with the uniform mag- 
netic field H is irrelevant. Thus, at T = the system 
remains critical for a certain range of the Ising field H. 
In the renormalization scenario outlined in Rcf. 10] a 
phase transition to the long-range ordered phase occurs 
when H grows large enough. In the SOS language this 
ordered phase is flat, and transition is of the 'roughening' 
type, by duality related to the Kosterhtz-Thouless phase 
transition p^ . 
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B. Renormalization flow 

For an analysis of the renormalization flow at nonzero 
magnetic field H it is necessary to include the Coulomb 
gas coupling because the field, although irrelevant at 
(/R = 2, tends to suppress height differences, i.e., to in- 
crease f/R. As deduced in Ref. 10] such an effect can 
also be realized without breaking the Ising symmetry, by 
the introduction of ferromagnetic next-nearest-neighbor 
interactions into the model Although we restrict our 
numerical investigation to the 2-parameter model |^ , the 
renormalization analysis still requires a set of 3 nonlinear 
scaling fields, which are chosen as 



t = 

h = 
ti = 



9 

„'2K 



(28) 



The constant ai determines the scale of t which remains 
to be determined. Apart from that, the expansion coef- 
ficients Uj are in principle unknown. The Ising temper- 
ature field ti is in lowest order chosen as the Boltzmann 
factor of an elementary Ising excitation, i.e., a triangle 
with three equal spins. It thus describes the activity of 
the magnetic charges m = ±1 in the Coulomb gas. In 
order to describe the renormalization of these parameters 
in the immediate vicinity of the fixed line h = ti ^ Q, the 
following renormalization exponents apply 



yt 

Vh 



9 
2^ 

2 



(29) 



First we address the special case T — ti — Q. Because of 
the marginality of t we add a nonlinear term in the flow 
for t. In differential form, the equations become 



Ah{l) 
d/ 



dt{l) 
"dT 



(30) 



where I parametrizes the rescaling factor h asb — exp(Z). 
In principle one has an unknown amplitude in the second 
equation, but we have disposed of it by a proper choice 
of ai. Thus, the scale of t is set such as to simplify the 
renormalization equations to (|30|l . The equation for h 
follows from the usual form h' = W'^h after substitution 
of b in terms of /, of yt using Eq. (|29|l . and 5r in terms 
of t. The sign in the equation for t follows because h 
suppresses the height differences in the SOS language. 
The flow equations H30|l are equivalent to those describing 
the Kosterlitz-Thouless and roughening transitions. 
Elimination of / from Eq. H30|l and integration yields the 
trajectory in the h, t plane as 



(31) 



where the constant c follows from the initial conditions 
which are chosen as /i(0) = 1/4-1- 5h, t{0) = 1/4. The 



physical motivation of this choice is that the KT transi- 
tion line obeys h = t, so that we select a point at a dis- 
tance 6h to the KT point of the nearest-neighbor model 
Eq. (Q. For small 6h one finds « 6h/2. Elimination 
of h in Eq. leads to 



di(0 _ _ ,2 „ 2 
dl ^ 



(32) 



Integration, substitution of the initial conditions and 
some rearrangement lead to the renormalization flow for 
ii = as 



(33) 



Next we introduce a nonzero Ising temperature T so 
that also ti > 0. It seems reasonable to assume that the 
renormalization flow of h and t is not seriously affected 
for small ti. For simplicity we make a stronger assump- 
tion, namely that the flow of h and t is independent of 
that of ti. We first focus on the question whether the 3- 
state Potts critical line in the {H, T) plane extends to the 
KT transition at T = 0, or, in other words, whether there 
are points in the immediate vicinity of {H — Hkt,T = 0) 
that flow towards a region where we can be confident that 
a Potts-like transition occurs. For small 6h/2 and ti the 
first part of the path will be governed by the KT fixed 
point. We assume that the flow will bring the system to 
a boundary which separates the regions governed by the 
KT and Potts fixed points. This boundary is obviously 
not determined in a quantitative sense, but this does not 
matter for the present scaling argument. Let it be sufh- 
cient to define this region by requiring that h and ti reach 
values of order 1 . 

We search for this region by choosing the (somewhat 
arbitrary) renormalized temperature = 3 where we 
have evidence |2^ that the Potts transition connects to 
the neighborhood oiT ^ H = 0. The shape of the critical 
line at ^r = 3 is determined by the flow equations for ti 
and h, namely 



h{b) = b^^Hi 
h{b) = b^^h 



(34) 



where one may take yti — yh — 1/2 as long as ti and 
h are small so that the change of g^i can be neglected. 
Then, t\{b)/h{b) is constant along a flow line and can be 
chosen such that the model is critical, say for 



ti{b)/h{b)^ 13. 



(35) 



where /? is a constant of order 1. For larger values of 
ti and h the relation will no longer be linear but it is 
reasonable to expect, and in agreement with numerical 
results [23. 125I I , that there is a fair range where /? is still 
of order 1 . 

Thus we consider a point /i(0) = 1/4 + <5/i, t(0) = 1/4, 
ti(0) = 8ti in the vicinity of the KT point and apply 
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FIG. 2: Sketch of the renormahzation flow in the parame- 
ter space of the renormahzed couphng gn, the scahng field h 
and the Ising temperature field ti. The flow of and h is 
anomalously slow near the point gn, = 9/4, /i = 0: most of 
the growth of ti occurs here. 



a transformation such that the system flows to gri = 3 
or t = —1/2. According to Eq. H33|l the scale factor 
of this transformation is b = exp{l) = exjp(Try^2/6h). 
It follows from Eq. (|33|l that (for small 5h and Sti) the 
system is located near the KT fixed point h — 0, t ^ Q 
for most of the range of I. Therefore, the flow of the 
Ising temperature field ti is determined by the exponent 
ytj = 7/8 at = 9/4 or t = 0. Thus, at = 3 it reaches 
the value tj = b'^^^Sti = exp(77r/4V25X)5ii. Since h = 
— t = 1/2 up to unimportant corrections, Eq. (|35|l leads 
to 



dti — — exp ; 

2 AV2ShJ 



(36) 



which solves dti for all Sh > 0. This implies that the 
Potts critical sheet connects to the KT point. The 
resulting renormahzation flow is sketched in Fig. |21 
Substitution of Eq. and H — Hkt + SH leads in 
lowest order to 



2 



2 4v/2a(iJ - Hkt) 



(37) 



where a is a function of the aj . This equation determines 
the shape of the critical line near the KT point of the 
model and is clearly incompatible with the prediction 
of the invariant-polynomial scenario, Eq. (|12|1 . 



axis determines the transfer direction. We have used two 
choices for the orientation of the lattice: one set of bonds 
parallel or perpendicular with respect to the axis. For 
the first case one may apply a decoration transformation 
to one half of the parallel bonds in order to construct a 
symmetric transfer matrix. However, the decoration of 
antiferromagnetic bonds leads to complex weights which 
we wish to avoid. We have thus used a non-symmetric 
transfer matrix, in combination with a suitable tridiago- 
nalization method to find the leading eigenvalues. These 
were obtained for even linear system sizes up to L = 22, 
which corresponds with an actual finite size of llV3 
nearest-neighbor bonds. The second construction, with 
a set of edges perpendicular to the transfer direction, 
leads to a symmetric matrix when two layers of spins 
are added. This allows the use of the conjugate-gradient 
method which is, in our applications, more stable than 
the tridiagonalization method. Finite-size calculations 
with L multiples of 3 up to L — 24 were performed using 
this second construction. 

A sparse-matrix decomposition was used for both con- 
structions. Most of the technique is already implicit in 
the work of Nightingale 0. Further details are listed 
in Ref. |126|| which concerns the case of the honeycomb 
lattice, but the essential steps are applicable to the tri- 
angular lattice as well. During the analysis of the re- 
sults of both types of transfer matrix we found that they 
were mutually consistent. Furthermore it became clear 
that the second transfer matrix, with a set of bonds per- 
pendicular to the transfer direction, allowed a somewhat 
more accurate analysis. In the following we describe the 
situation of the second construction. 

For T = the transfer matrix decomposes in a num- 
ber of diagonal submatrices characterized by a conserved 
number of 'strings' so that the numerical diagonalization 
task simplifies. The transfer-matrix construction for this 
case has been outlined in Ref. and enabled the study 
of systems with linear sizes up to L = 27. 

The magnetic correlation function along the coordinate 
r in the length direction of the cylinder is defined as 
gmif) = (soSr)- At large r, this correlation function 
decays exponentially with a characteristic length scale ^ 
that depends on K, H and L 



g„iir) cx e 



^-r/aK,H.L) 



(38) 



and can be calculated from the largest two eigenvalues 
Aq and Ai of the transfer matrix: 



IV. NUMERICAL METHODS 



(39) 



A. Transfer-matrix calculations 

Most of the the transfer-matrix calculations were per- 
formed for T > so that we had to use a binary repre- 
sentation for the Ising spins, leading to a transfer matrix 
of size 2^ X 2^ for a system with finite size L. We de- 
fine the spin lattice on the surface of a cylinder, whose 



where the factor \/3 is a geometric factor, i.e., the ratio 
between the thickness of two layers added by the transfer 
matrix and the length of a nearest-neighbor bond. The 
significance of these relations lies in the fact that the 
assumption of confornial invariance links ^ on the 
cylinder with the magnetic scaling dimension (one 
half of the magnetic correlation function exponent rj). In 
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terms of the scaled gap 



X„,{K,H,L) 



L 



(40) 



one has Xm{K, H, L) ~ Xm in the hmit of large L. Since 
the three-state Potts universal value of the magnetic scal- 
ing dimension is known to be Xm = , and the transfer- 
matrix algorithm evaluates Xm as a function of its argu- 
ments, one can find a numerical approximation to the 
critical value of K for a given value of H or vice versa. 
The shape of the critical line prescribes the use of differ- 
ent ways in different regions. For small H and large \K\, 
the critical line is almost parallel to the zero-field line, so 
that it becomes more efficient to solve for H than for K. 

As a consequence of corrections to scaling, the solu- 
tion will not precisely coincide with the critical point. 
The effects of an irrelevant scaling field u and a small 
deviation t with respect to the critical value oi H oy K 
are expressed by 



Xm{K, H, L) = Xm + auD>^ -f 



(41) 



where a and h are unknown constants, Xm — 2/i — 
and y* = I for the 3-state Potts universality class. Thus 
the solution for K of 



Xm{K,H,L) = 



15 



(42) 



which is denoted Kc^^\h, L), depends on the finite size 
L and the irrelevant field as 



Kc'^^\H,L)=Kc + ciLy^-y' + ■■■ 



(43) 



because the two correction terms in Eq. (|41|l must cancel 
and t (X Kc^^\H, L) — Kc- We thus generated sequences 
of iterated estimates of Kc by solving kJ'^\h, L) and 
ci{L) in the equations 

kJ^Hh,L) ^ Kc^'\H,l) + c,iL)iy'~y' (44) 

ioi I — L and L + 1. These sequences appear to converge 
faster with increasing L than the Kc'-^\h,L). Remain- 
ing corrections may be due to additional contributions to 
Eq. H41I) . for instance scaling as L^^'. Thus we defined 
Kc^{H,L) by solving the equation 



(45) 



for / — L and L ~\- \. Further estimates can be obtained 
with correction exponents 2j/i — 2yt, or by treating the 
correction exponents as a free variable in which case three 
values of I have to be used. Several variations of this 
procedure were tried which leads to some insight in the 
numerical inaccuracies of the fitting procedure. 

Our final estimates of the critical points for T > are 
listed in Tab. The apparent accuracy of the critical 
points is satisfactory for most of the field range, but it 
deteriorates rapidly at small fields. Nevertheless this re- 
gion has our special interest: we wish to determine how 



the Potts line connects to the KT point on the T = line, 
because this is where the theoretical predictions (see Sec- 
tions ^1 and ^Oj are markedly different. 

We have also reconsidered the determination of the 
KT point at T = given in Ref . (Ij] . In that work, the 
finite-size data for the critical field i?KT were obtained 
by requiring that the scaled gap associated with the spin 
waves of period 6, i.e., electric charges e = ±1 were equal 
to the expected value Xi_q = 2/9 at the KT transition. 
These estimates HxriL) of the KT point, which were 
obtained for system sizes up to L = 27, were found to be 
considerably size-dependent. They were fitted accord- 
ing to i?KT(i) = Hkt -I- a/(& -I- InL) -I- cL~^ which led 
to extrapolated estimates Hkt that displayed only a re- 
markably small size-dependence. On this basis, the final 
estimate was given as Hkt = 0.266 ± 0.010 in Ref. [T^ . 
In our present work we have reproduced these data for 
Hkt{L). We have not used the procedure that requires 
that the scaled gap is equal for two subsequent system 
sizes. Since all points in the range H < Hkt satisfy 
this scaling equation asymptotically, the solutions may 
converge to any point in this range, depending on the 
corrections to scaling. This is another reason behind the 
discrepancy mentioned in Sec. ^concerning an earlier re- 
sult for the location of the KT point ■ We have used 
the Hkt{L) data as input for several other iterated fit 
procedures consisting of subsequent extrapolation steps 
according to behavior, and powers of 1/lnL. These 
fits led to results for Hkt that were rather consistently 
close to 0.26, with differences up to 0.02. 

As an independent approach we have estimated Hkt 
from the requirement that the scaled gap based on the 
magnetic dimension Xq i is equal to the expected value 
1/8 at 5r — 9/4. Since such fractional magnetic charges 
(corresponding with vortices of strength 2) do not exist in 
this model, this scaled gap cannot directly be calculated 
for a fixed system size L. However, it can be obtained by 
combining free-energy data for system sizes L = 3n ± 1, 
where n is an integer, as explained in Ref. jl^ . The 
same extrapolation procedures as above were tried, and 
led to results again consistent with Hkt = 0.26, but 
with differences up to about 0.04. Our final conclusion 
is Hkt = 0.26 ± 0.02, similar to the value presented in 
Ref. |l2l | but with a slightly more conservative error esti- 
mate. 

The numerical results for the critical points are com- 
bined in Fig.m the phase diagram in the {H, e^^) plane. 
We remark that for large \K\ and relatively small H, the 
solutions for H become strongly finite-size dependent and 
slowly convergent. This problem is apparently due to the 
proximity of the KT transition at K — ~oo. This is il- 
lustrated in Fig. 13 in which a set of lines represent the 
finite-size solutions for L = 3, 6, • • • , 24 together with 
the extrapolated critical line. On the basis of our limited 
range of finite system sizes, the estimation of the critical 
points thus becomes increasingly difficult for large \K\. 

In order to provide further justification for our assump- 
tion that the critical line belongs to the 3-state Potts 
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TABLE I: Extrapolated results for selected points on the 
critical line. 



# 


H 


K 


# 


H 


K 


1 


0.55(5) 


-2.3 


20 


1.85 


-0.759438 (2) 


2 


0.57(3) 


-2.2 


21 


1.90 


-0.755049 (1) 


3 


0.59(2) 


-2.1 


22 


1.95 


-0.751498 (2) 


4 


0.610(10) 


-2.0 


23 


2.0 


-0.748715 (2) 


5 


0.634(4) 


-1.9 


24 


2.1 


-0.745199 (1) 


6 


0.658(2) 


-1.8 


25 


2.178 


-0.744130 (1) 


7 


0.6885(10) 


-1.7 


26 


2.3 


-0.744958 (1) 


8 


0.7219(3) 


-1.6 


27 


2.4 


-0.747586 (1) 


9 


0.7607(1) 


-1.5 


28 


2.5 


-0.751708 (1) 


10 


0.8 


-1.414 (1) 


29 


2.75 


-0.7673233(4) 


11 


0.9 


-1.2395 (2) 


30 


3.0 


-0.7887774(4) 


12 


1.0 


-1.11422 (1) 


31 


3.25 


-0.8145143(3) 


13 


1.1 


-1.02100 (1) 


32 


3.5 


-0.8434661(2) 


14 


1.2 


-0.95030 (1) 


33 


4.0 


-0.9082113(2) 


15 


1.3 


-0.896040 (5) 


34 


4.5 


-0.9791030(2) 


16 


1.4 


-0.854175 (2) 


35 


5.0 


-1.0539340(2) 


17 


1.5 


-0.821890 (1) 


36 


5.5 


-1.1313743(2) 


18 


1.6 


-0.797164 (2) 


37 


6.0 


-1.2105830(3) 


19 


1.75 


-0.771064 (3) 










FIG. 3: Finite-size solutions for the critical points, and our 
final estimates in the region of small field and temperature. 
The dashed lines connect the solutions shown as triangles (A). 
From right to left the lines show data for finite sizes L = 3, 
6, ■ ■ ■ , 24. The solid line with circles (o) indicates our final 
estimated result for the critical line. 



universality class (used for the determination of the crit- 
ical points) we perform a consistency test by calculating 
the conformal anomaly c at the estimated critical points 
at H = 1.5, 2.5, 3.5 and 4.5. Iterated fits similar to 
those used for the calculation of the critical points were 
applied. All these results are consistent with the exact 
value c = I . The error margin varies between a few times 
10"'^ for H = 1.5 and a few times 10~^ for H — 3.5. In 
comparison with previous work these results further 
restrict the scale of possible deviations from 3-state Potts 



universality. 

B. Monte Carlo results 

The apparent difficulty to obtain accurate critical 
points for small H by the transfer-matrix method invites 
further investigation by means of Monte Carlo simula- 
tions, which allow the use of much larger finite sizes. 
In particular we determine how the critical amplitudes 
behave for small T, and make a comparison with the 
renormalization prediction. To this purpose we define 
a specific- heat-like quantity C, i.e., the second order 
derivative of the free energy to a parameter conjugate 
to an energy-like density in the Hamiltonian, for which 
we may take the magnetization. Indeed the Ising field H 
drives the Potts-like transition to the ordered state (ex- 
cept at the maximum of the T vs. H curve) and thus 
plays the role of the temperature in the Potts model. We 
thus define C by 

C=^-N{{m')^(mf) (46) 

where m is the Ising magnetization. Similarly, we define 
a quantity similar to the magnetic susceptibility of the 
Potts model. In the q = 3 Potts model, the zero-field 
magnetic susceptibility can be expressed in magnetiza- 
tion fluctuations by 

X = N{ml) (47) 

where rup — n\ + + n| — nin-i — n2n^ — ny,ni expresses 
the Potts magnetization in terms of the densities ni of 
Potts variables in state i. In the scaling context of the 
present scaling analysis, the densities rn may be defined 
as the number of minus-spins on sublattice i. Thus x 
describes the response of the model to staggered fields 
acting on the Ising spins. 

The simulations used triangular L x L lattices with 
periodic boundary conditions. We used a combination 
of the standard Metropolis algorithm and the geometric 
cluster method 28]. The latter method executes nonlo- 
cal updates and leads to a faster relaxation. But it does 
not change the Ising magnetization. For this reason also 
Metropolis steps were included. First we sampled C in a 
suitable range of H, to study its divergence, at fixed val- 
ues of coupling K which are taken from Table The re- 
sults are shown in Figs. 0] and [Hand show that the finite- 
size divergence of C at the critical line becomes weaker 
when H decreases. In order to study this phenomenon 
in a more quantitative sense, we have determined C and 
X at several critical points taken from Tabled for several 
system sizes. Results for H — 0.61, 0.658, 0.8, 1.0, 1.5, 
2.0 and 5.0 are shown in Fig. Eland Fig. [7| 

Finite-size scaling of the free energy density of a system 
with finite size L can be expressed by 

f{t,h,u,L) = L-'^f{LyH,Ly'^h,Ly^u,l)+g{t,h,u) (48) 
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0.35 



0.3 



U 0.25 



0.2 



0.15 




0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 
H 

FIG. 4: Specific heat-like quantity C versus field H for two 
values of the Ising coupling, K — —1.0 (o) and —1.5 (□). 
Both sets of lines display data for system sizes L = 12, 24, 
48, and 96. These two cases indicate that the amplitude of 
the divergence of C decreases with decreasing field H. 

0.8 I ^ r- ^ , , 1 




0.1 ' ' ' ' ' ' 1 

4.8 4.9 5 5.1 5.2 5.3 5.4 

H 

FIG. 5: Specific heat-like quantity C versus H at Ising cou- 
pling K = —1.053934, near the critical point &t H = 5.0. 
In comparison with Fig. |1| the finite-size divergence is much 
stronger. 

where i, /i, and u denote the temperature, magnetic field 
and irrelevant field respectively, and g the regular part of 
the transformation. Differentiation of / yields the scaling 
behavior of the quantities C and x ^ 

C{u, L) = Co+ L^y'-'^c{Ly^u, 1) 



= Co + L^y*-\bo + biL'^u + h2L^y'u^ + ■■■) (49) 

and 



= Xo + L^y^-^ibo + biLy^u + b^L'^y^u^ + ■■■) (50) 

We have fitted the numerical data by these two equa- 
tions and thus derived the amplitudes listed in Tab. ^ 



0.2 




2 g L f f f 

1.5 2.5 3.5 4.5 5.5 

InL 



FIG. 6; Specific heat-like quantity C versus system size on 
a double logarithmic scale. Data are shown for system sizes 
6 < L < 192 at seven points [H, K) on the critical line. The 
symbol • represents H = 0.61; V: H = 0.658; x: H = 0.8; 
A: H = 1.0; n: H ^ 1.5; O: H ^ 2.0; m: H ^ 5.0. The seven 
lines represent the fitted results. The statistical errors are not 
shown in this figure. They do not exceed the thickness of the 
lines. 
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2 ' ' ' ' 1 

1.5 2.5 3.5 4.5 5.5 

InL 

FIG. 7: Susceptibility-like quantity x versus system size L(= 
6, • • ■ , 192) on a double logarithmic scale. Data are shown 
for seven points {H, K) on the critical line. From bottom to 
top: H = 0.61, 0.658, 0.8, 1.0, 1.5, 2.0, and 5.0 respectively. 
The two lines for H = 2.0 and 5.0 coincide on this scale. The 
statistical errors are not shown in this figure. They do not 
exceed the thickness of the lines. 

and Tab, mil The amplitude b^ of the leading divergence 
of C decreases with field except close to the maximum 
of the critical line in the (7J, T) plane. At the maxi- 
mum the field H fails to bring the system into the or- 
dered phase and the amplitude 60 thus vanishes. Also 
for the susceptibility-like quantity x the amplitude &o of 
the finite-size divergence, shown in Tab. IIIII decreases 
regularly when the KT point is approached. 

The behavior of the amplitude 69 at small field and low 
temperature follows from the renormalization-flow analy- 
sis. Starting from a point in the vicinity of the KT point, 
we renormalize until we arrive at the boundary with the 
region dominated by the 3-state Potts fixed point. Let 
&KT be the corresponding scale factor. Since the specific 
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TABLE II: Parameters describing the finite-size behavior of 
the specific heat-like quantity C. The third column is the 
amplitude 60 of the leading divergent term. The amplitudes 
for H = 1.5, 2.0 are relatively small because the critical line 
runs almost parallel to the field direction. The column cor- 
responding with bi {i = 1,2) are the irrelevant corrections 
amplitudes. 



H 


Co 


bo 


bi 


b2 


0.61 


0.188(2) 


0.0026(2) 


0.027(3) 




0.658 


0.196(2) 


0.0039(3) 


0.007(3) 




0.8 


0.192(3) 


0.0126(4) 


-0.022(4) 




1.0 


0.151(3) 


0.0285(6) 


-0.035(5) 




1.5 


0.097(2) 


0.0254(4) 


-0.055(3) 




2.0 


0.084(3) 


0.0029(3) 


-0.015(6) 


-0.097(2) 


5.0 


-0.096(6) 


0.131(2) 


-0.031(9) 





TABLE III: Parameters describing the finite-size behavior of 
the susceptibility-like quantity x- The third column is the 
amplitude bo of the leading divergent term. It decreases reg- 
ularly with the field. The next columns shows the irrelevant 
correction amplitudes bi (i — 1,2). 



H 


Xo 


bo 


bi 


ba 


0.61 


-17.5(18) 


0.5281(12) 


-0.143(50) 


16.3(16) 


0.658 


-10.5(20) 


0.5896(14) 


-0.142(55) 


2.88(36) 


0.8 


-18.7(37) 


0.7675(24) 


-0.75(10) 


18.1(33) 


1.0 


0.95(11) 


0.9046(17) 


-0.164(25) 




1.5 


0.41(10) 


1.1209(12) 


-0.103(24) 




2.0 


0.10(14) 


1.2130(24) 


-0.003(34) 




5.0 


-0.51(14) 


1.2133(22) 


0.070(37) 





heat-like quantity C is defined by means of differentia- 
tion of tiie free energy to tiie uniform field H , we keep 
track of how H changes under this transformation. The 
marginality of Sh at = 9/4 is expressed by Eq. (|31|l : 
when we write h — t + 6h, it is clear that Sh varies only 
by a factor of order 1 as long as t is of order 1 . In the con- 
text of scaling, we thus have Sh' « Sh where the prime 
indicates the value at the boundary. Within the Potts 
region we rescale the system to size 1 with the remaining 
scale factor L/6kt 



Sh" 



(— 

\bKT 



Sh' 



(51) 



where the Potts temperature exponent | applies because 
it corresponds with the Ising magnetic field. The behav- 
ior of C follows as 



h,\)^Ln^^ r 



const 



5V 

(52) 

The susceptibility-like quantity x is obtained by dif- 
ferentiation of / to the staggered Ising field, which, as 
explained in Ref . [Tol | , is associated with a Gaussian spin 



wave perturbation of period p = 6, i.e., with electric 
charges e = ±1. The exponent of the staggered field 
thus takes the value 2 — Xi q = 16/9 at the KT fixed 
point. Therefore, at the boundary with the Potts region 
we have 



/ 16 



(53) 



Within the region dominated by the Potts fixed point, the 
magnetic exponent = 28/15 applies. Renormalization 
with the remaining scale factor L/6kt leads to 



L 
6kt 



(54) 



so that X scales as 



dhl. 



dhl 



f{LTSh^fh,,,l) = 

2fi 8 

LT-.blff" 



const 



(55) 



According to Sec. IIIIBI rescaling by a factor 6kt re- 
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suits in an Ising temperature field b^^^e 



= /3/2 so that 

^KT oc e^~^ . For strong coupling the renormalization 
scale 6kt is large, which is indicative of the crossover 
phenomenon close to the KT transition. There we need 
large system sizes L > 6kt in order to reach the vicinity 
of the Potts fixed point. The substitution of 6kt into 
Eq. llll) and Eq. leads to 



C = Lte^^ f" + const 



and 



T — 



26 128 ,, 
— 3T5« J" 



const 



(56) 



(57) 



From a comparison of Eqs. (|49|l and 1)56(1 . and of Eqs. 
and (jSTJ, we expect that 5o oc e^-^ for C, and 60 (x 
e3i5^ for X, when \K\ is large enough. We thus expect a 
linear relation between In 60 and K for sufficiently strong 
coupling K. A fit to the numerical data yields the slopes 
as about 2.8(1) for C and 0.64(4) for x- These spoles 
do not agree accurately with the analytic values. This 
suggests that the Ising temperature used for the calcula- 
tion of the amplitudes is not small enough. However, the 
qualitative amplitude dependence is reproduced, and the 
rough agreement suggests that we are not far away from 
the asymptotic regime. 



V. FIT AND DISCUSSION 

A. Roots of the 0(2) invariant polynomials 

The number of 37 critical points in Table |l] together 
with the additional KT point (/i, z) = (/ikt!~1) is suf- 
ficient to attempt fits of 0(2) invariant polynomials 1191 
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.1 I . . . i 

0.2 0.4 0.6 0.8 1 

tanh H 

FIG. 8: Phase diagram in the {h, z) plane (where h — tanh JJ" 
and z = tanhA"). The KT point is denoted as ■, and the LG 
point, which is Baxter's hard-hexagon model, as ». 



up to order 10. We have not used data points 1 and 2 
because of their limited accuracy, and performed a least- 
squares fit to the remaining 35 data points for z > — 1. 
We have also tried direct fits to several subsets of these. 
In each case we have investigated the effect of enforcing 
the curve to pass through the KT point, or to extract 
the value of Hkt from the fit. It is found that the co- 
efficients in the equation — are not flexible enough 
to even qualitatively fit the numerical data. The least- 
squares fit to /g = excluding the KT point consists of 
two avoiding solutions, which lead to 2 disconnected 'crit- 
ical' lines which have unphysical ranges. For z — — 1 one 
line terminates at h* « 0.501 > Hkt- When enforcing 
the KT point at /ikt = 0.25, the two avoiding branches 
repel one another even stronger. Direct fits to different 
subsets of critical points lead to similar results. Fits to 
higher order equations fs — and /lo = display the 
same problems. Even more avoiding solutions enter. The 
numerical problems are clearly displayed by the values of 
the fitted coefficients, which span a range of many orders 
of magnitude. In summary, the roots of invariant poly- 
nomial equations cannot fit the critical curve. The main 
problem is the approach of the curve to the KT limit im- 
posed by Eq. H12|l : all such roots approach the KT point 
vertically in the {h, z) plane, whereas the numerical data 
in Fig. |S1 indicate a horizontal approach. 

Sufficiently far away from the KT point, the problem 
dissolves, and our numerical data for the critical points 
can well be approximated by means of invariant polyno- 
mials. For example, the polynomial of order 6 can re- 
produce the critical points for H > 1.5 within the error 
margins quoted in Table U The coefficients, determined 
by means of a least-squares fit, are listed in Tab. II VI 



B. The renormalization solution for small field 

For small field we expand Eq. (|37|) and take into ac- 
count higher order terms in the physical fields. This leads 



TABLE IV; Coefficients {i = 0, ...,6;j = 1,2,3) of the 
invariant polynomial fa, Eq. 11911 . The condition for criticality 
reads — 0. 







^^2 


-0.002942307242 


-0.097656582607 


2.725303204552 






4 


-0.000000000017 


0.000351055639 


-0.020534345140 




2 




0.044511873969 


0.660086779714 


0.000001077354 






4 


-0.002865984972 


-0.084578467922 


1.000000000000 


1.4 1 . . . ^ 



1.2 - 



0.8 - 

' 0.6 - jT 

0.4 - /? 

0.2 - / 

I i . . . ^ 

0.5 1 1.5 2 

H 

FIG. 9: Phase diagram in the (H,K~^) plane. The symbol 
o denotes the Potts transition points; ■ the KT point; and 
the solid line describes the fit of the expression based on the 
renormalization prediction for the critical line in the small 
field region. 



to 

-1= E «.(^-^KtF^' (58) 

The numerical data for the critical points for H < 1.75 
are fitted satisfactorily {i.e., within the error margins 
quoted in Table P) by this formula using 6 coefficients. 
The numerical results and the fitted function are shown 
in Fig. 1^1 The values of the coefficients are listed in 
Tab. El 

TABLE V: Coefficients aj of the expansion 15811 of ~K^^ in 
powers of [H — Hkt)^ ■ 



ai 


a2 


as 


0.487432 


0.119116 


0.765066 


a4 


as 


0.6 


1.017104 


-1.949253 


0.652161 
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VI. CONCLUSION 

The invariant-polynomial scenario formulated in 
Sec. m and the renormalization scenario formulated in 
Sec. mil lead to analytic expressions for the critical line 
in the (iJ, T) diagram that are mutually inconsistent for 
T I at finite H . This shows that at least one of the 
underlying assumptions must be incorrect. The renor- 
malization prediction appears to successfully describe the 
numerical data for small \H\. Although the asymptotic 
regime is not quite reached (as can, for instance, be seen 
in Fig. ini where the leftmost points behave almost lin- 
early instead of as a square root), an asymptotic expan- 
sion leads to an accurate description of the data, and 
allows a smooth extrapolation to zero Ising temperature 
in agreement with Eq. (|37|l . The analysis in Sec. IIV Bl 
of the critical amplitudes fits precisely in this picture. 
Thus our analysis does not give reasons to doubt that 
the renormalization scenario correctly describes the es- 
sential physics of the model near the KT transition. 

In contrast, the invariant-polynomial scenario does not 
agree with the numerical data. It predicts a 'vertical' ap- 
proach to the KT point in the (/i, z) diagram (see Fig.|SI) 
where it should be horizontal. Our interpretation is that 
the assumption of analyticity of the critical line in the 
(/i, z) parametrization is false at the KT point, so that 
the line cannot be described by the zeroes of a polynomial 



of a finite order. 

Since it now appears that the invariant-polynomial sce- 
nario fails in the case of the triangular-lattice Ising model, 
the question arises whether similar, apparently success- 
ful, analyses of the critical lines of the honeycomb- and 
the square-lattice Ising model in terms of invariant poly- 
nomials [T^ I20I I have to be reconsidered. Here we may 
point at the simpler topology of the (iJ, T) diagram for 
the honeycomb and the square lattices: the critical line 
connects to T = only in the lattice-gas points H = ±00. 
In the case of the triangular lattice model, crossover phe- 
nomena near the KT point are responsible for the nonan- 
alytic 'shape' of the critical line. In the absence of such 
crossover phenomena, there is no inconsistency with the 
invariant-polynomial scenario, and our present analysis 
has therefore no direct consequences for the work pre- 
sented in Refs. [Sll^. 
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